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ABSTRACT 

Structural parameters (half-light radius r e , mean effective surface brightness 
(/i) e , and Sersic index n, parameterizing the light profile shape) are derived for a 
sample of galaxies in the rich cluster AC 118 at z — 0.31: so far the largest (N=93) 
sample of galaxies at intermediate-redshift with structural parameters measured 
in the near-infrared. The parameters are obtained in two optical wavebands 
(V and R) and in the K —band, corresponding approximately to B, V and H 
rest-frame. 

The distributions of r e at z — 0.31 match those for the Coma cluster (i.e. for 
the local universe) both in the optical and in the NIR. The K —band distribution 
is of particular interest, since the NIR light mimics the mass distribution of 
galaxies. The similarity of the distributions for the two clusters (AC 118 and 
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Coma) proves that the galaxies at the bright end of the luminosity function did 
not significantly change their sizes since z ~ 0.3 to the present epoch. 

The ratio of the optical to the NIR half-light radius shows a marked trend 
with the shape of the light profile (Sersic index n). In galaxies with n > 4 (typical 
bright ellipticals) t^nir ~ 0 . 6 r e , O pt , while the average ratio is 0.8 for galaxies with 
lower n (typical disk systems). Moreover, the NIR Sersic index is systematically 
larger than in the optical for n < 4. These results, translated into optical and 
optical-NIR color gradients, imply that the optical color gradients at z ~0.3 are 
similar to those of nearby galaxies. The optical-NIR color gradients are in the 
average larger, ranging from -0.73 mag/dex for n < 4 to -0.35 mag/dex for n > 4. 

Models with ‘pure age’ or ‘pure metallicity’ gradients are unable to reconcile 
our color gradients estimates with observations at z ~ 0, but we argue that the 
combined effects of age and metallicity might explain consistently the observed 
data: passive evolution (plus the possible effect of dust absorption) may account 
for the differences between the optical and NIR structural properties. 

The lack of any major change in r e ,Nm since z ~ 0.3 suggests that merging 
involving bright galaxies did not play a significant role in the last rN " , 4.4 Gyr 
(Om= 0.3, = 0, H 0 = 50 km s _1 Mpc _1 ). 

The results of the present paper will be applied to the study of the scaling 
laws in subsequent works. 

Subject headings: galaxies: clusters: individual: AC 118 - galaxies: evolution 
- galaxies: fundamental parameters (effective radius, surface brightness, Sersic 
index ) - galaxies: photometry 


1. Introduction 

Our knowledge on the evolution of galaxies comes from the study of the properties of 
families of galaxies at different redshifts. Observable quantities like magnitudes, colors, struc¬ 
tural parameters, velocity dispersions and line strengths are combined into well-established 
relations like, for instance, the color-magnitude relation (e.g. Bower, Lucey, & Ellis 1992; 
Kodama et al. 1998), the Mg-cr relation (Ziegler & Bender 1997), the fundamental plane 
(Dressier et al. 1987; Djorgovski & Davies 1987; Jprgensen, Franx, & Kjaergaard 1996) and 


1 Based on observations collected at European Southern Observatory (ESO n. 62.0-0369, 63.0-0257, 
64.0-0236) and on data from the STScI Science Archive 
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its projection in the r e - (/x) e plane, that is the Kormendy relation (Kormendy 1977; Capac- 
cioli, Caon, D’Onofrio 1992; Ziegler et al. 1999). These relations can now be studied over a 
significant fraction of the age of the universe. Studies of the properties of early-type galax¬ 
ies in rich clusters lead to a large amount of observational evidence in favor of the passive 
evolution scenario, in which an early (z >2) and short (~1 Gyr) burst of star formation is 
followed by a long period where the aging of stellar populations drives the observed proper¬ 
ties of galaxies through intermediate redshifts until the present epoch (Kodarna & Arimoto 
1997; Stanford, Eisenhardt, & Dickinson 1998; Jprgensen et al. 1999; Kelson et al. 2000a, 
and references therein). On the other side, the alternative scenario of the hierarchical merg¬ 
ing seems also consistent with most of the existing data (e.g. Kauffmann & Chariot 1998), 
indicating either that the present observational evidence is not sufficient to discriminate 
between the two pictures or that galaxy evolution is driven by both mechanisms. 

The half-light radius and the mean surface brightness of a galaxy measure the size and 
the density of the luminous matter. These quantities are related both to the gravitational 
potential of the galaxy and to the properties of the stellar populations. The relations existing 
between the structural parameters of galaxies have been recently studied at intermediate 
redshifts, both in the form of r e vs. (/i) e (Ziegler et al. 1999, and references therein) and as 
luminosity-size relation (Schade et al. 1996; Simard et al. 1999). Schade et al. (1996) found 
that disks at z >0.5 are in the average ~ 1.5 mag brighter than in local galaxies. This was 
contrasted by the results of Simard et al. (1999) who found no luminosity evolution in disks 
up to z ~0.9. Ziegler et al. (1999) found that the Kormendy relation of cluster galaxies 
at z = 0.4 and z = 0.55 is consistent with the passive evolution scenario. The Kormendy 
relation has also proven to be a powerful tool to test the expansion of the universe by the 
measurement of the cosmological dimming (i.e. for the Tolrnan test). Lubin & Sandage 
(2001) studied the Kormendy relation up to z — 0.92, finding that the data are consistent 
with the cosmic expansion together with passive evolution. Finally, the structural parameters 
enter the fundamental plane of early-type galaxies, of which the Kormendy relation is a two- 
dimensional projection. Due to its very small scatter, the fundamental plane is a powerful 
tool to study the evolution of the M/L ratio of galaxies belonging to distant clusters. The 
recent determinations of the fundamental plane at intermediate redshifts constrain the epoch 
of star formation in early-type galaxies at large redshifts (z >2 for Q M —0.3 or z >5 for fi M =l 
and Da = 0; see van Dokkum & Franx 1996; Jprgensen et al. 1999; Kelson et al. 2000a). 
However, these studies have to face with the problem of the progenitor bias, i.e. with our 
ignorance of the relationship between local galaxies and their high redshift counterparts (see 
van Dokkum & Franx 2001 for a thorough discussion). 

A clue for the interpretation and for the applications of the scaling laws is the depen¬ 
dence of the structural parameters on wavelength. The basic distinction between optical and 



NIR is that the optical wavebands are very sensitive to the effects of metallicity (through 
line-blanketing) and to the content of younger populations (dust absorption included), while 
the NIR wavebands are dominated by the old, more quiescent, stellar populations and are 
virtually unaffected by dust absorption. For these reasons the NIR is the best tracer of 
the mass distribution inside a galaxy. Local galaxies are smaller at NIR wavebands than in 
the visible, with a similar, though weaker, trend also within the optical wavelength range 
(Scodeggio et al. 1998; Pahre, de Carvalho, & Djorgovski 1998). This behavior is con¬ 
firmed by the relatively large optical-NIR color gradients measured in early-type galaxies, 
the optical-optical gradients being usually much smaller, except when the U— band is in¬ 
volved (Peletier et al. 1990a; Pcletier, Valentijn, & Jameson 1990b; Hinkley & Im 2001; 
Nelson et al. 2001). Small optical color gradients were recently measured in distant (up to 
z ~1) cluster early-type galaxies by Saglia et al. (2000) and by Tamura et al. (2000). These 
authors concluded that metallicity is the primary factor driving the color gradients. 

The derivation of the structural parameters at intermediate redshifts is affected by the 
problem of the small angular size of the galaxies, whose half-light radii are typically less 
than 1 arcsec. For this reason many works were based on imaging taken with the HST. It 
has been demonstrated, however, that reliable structural parameters can be derived from 
ground-based photometry as well, on high resolution (< 0.1"/pxl) and exceptional seeing 
(< 0.4") conditions, provided that the properties of the PSF are carefully taken into account 
(Schade et al. 1996; Jprgensen et al. 1999). 

In this paper we obtain the first multi-band large set of structural parameters for galaxies 
at an intermediate redshift (z = 0.31). The data are derived from ground-based (ESO-NTT) 
photometry in two optical (V and R) and one NIR (K) wavebands. The derived values for 
the structural parameters will be published in a forthcoming paper addressed to the study 
of galaxy evolution by means of the scaling laws (Busarello et al. 2001a, in preparation). 

The paper is organized as follows. The sample and the data used in this study are 
presented in § 2. The derivation of the structural parameters is outlined in § 3. The 
structural parameters from HST imaging are derived in § 4, where the results from Sersic 
models and bulge-disk decomposition are compared. In § 5 we use HST photometry for 
a subset of galaxies to verify the reliability of ground-based structural parameters. The 
structural parameters are presented in § 6, where we make a comparison between optical 
and NIR properties at z ~ 0.3 and with the values at z ~ 0. In § 7 we analyze the optical 
and NIR structure of galaxies in terms of internal color gradients. The results are discussed 
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2. Description of the Sample 

AC 118 (Abell 2744) is a rich (‘Coma-like’) galaxy cluster at redshift z = 0.31 that was 
studied by several authors in relation to its galaxy population, X-ray and radio properties and 
gravitational lensing (Barger et al. 1996; Couch et al. 1998; Stanford et ah 1998; Williams, 
Navarro, & Bartelmann 1999; Wn, Xne, & Fang 1999; Allen 2000; Govoni et ah 2001). 

New photometric data and medium resolution spectroscopy has been collected during 
three observing runs (October 1998 - September 2000) at the ESO New Technology Telescope. 
The data of interest for the present work are the imaging in the V and R wavebands obtained 
with the EMMI instrument and the imaging in the K waveband obtained with SOFI. The 
photometry is complemented with HST archive data consisting of a single pointing on the 
cluster center. The relevant information on the data is given in Table 1. The reduction of 
the K —band image is described in Andreon (2001), while the reduction of the NTT optical 
data will be presented elsewhere (Busarcllo et al. 2002b, in preparation). The HST data 
were retrieved from the STScI Archive and reduced following the steps described in Couch 
et al. (1998). The NTT photometry was used to calibrate the HST imaging into the Johnson 
R —band as outlined in Appendix A. 

The samples used in the present study consist in: A) 50 galaxies in the HST field; B) 37 
galaxies in the HST field in common with the ground-based R —band sample; C) 58 galaxies 
in common in the V— and R— band; D) 93 galaxies in common in the R— and K— band. 

Samples C and D are the results of two selection steps. The first selection was based 
on photometric redshifts (Massarotti, Iovino, & Buzzoni 2001), as derived from the above 
photometry plus J— and U —band (new NTT data and NTT archive data respectively). The 
reliability of the photometric redshifts was verified on 75 galaxies with available spectra 
(Merluzzi et al. 2002, in preparation). Morphology is not a selection criterion in this study, 
since a morphological classification is not possible from our data. The photometric redshifts 
algorithm flags the galaxies according to whether their colors are consistent with early- 
or late-type template spectra. We did not use this classification as a selection criterion, 
and retained in the samples all the galaxies with the correct photometric redshift. We 
remark, however, that only 7 galaxies out of 93 in sample D are ‘late-type’ according to their 
photometric properties (from U to K). The second selection step is based on the reliability 
of the measurement of the structural parameters and is the result of the analysis presented 
in § 3.3. This lead to include in the sample the galaxies brighter than R =21 and H=20.5, 
and to exclude objects with an estimated r e smaller than 0.15 arcsec. 

Samples A and B were selected in order to obtain 1) reliable structural parameters in 
the NTT A—band and 2) the largest number of galaxies in the HST field with the same 
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range of sizes and luminosities of the ground-based samples. This lead to include galaxies 
brighter than R =21 and to exclude objects with r e smaller than 0.15 arcsec both in the NTT 
and in the HST image. A further (a posteriori) selection was done for all the samples on the 
basis of the results of the fitting: objects with very large uncertainties (>100%) were not 
considered in the analysis. This however lead to exclude only 2% of the objects. 

Morphological classification from HST imaging by Couch et al. (1998) is available for 21 
out of 93 galaxies in sample D (17 out of 58 in sample C). Out of these, 6 (6) are classified as 
ellipticals, 10 (7) as SO, and 5 (4) as spirals. These numbers reflect the overall distribution 
of morphological types in AC 118, which is dominated by early-type (mainly SO) galaxies 
(Couch et al. 1998). The question of morphological types will be resumed in the conclusions. 
The relevant information on the samples is summarized in Table 2. 


3. Derivation of the structural parameters 

Since the typical galaxy sizes at z ~0.3 are smaller than ~ 1", it is crucial to take into 
account the effects of the Point Spread Function (PSF) when deriving structural parameters 
from ground-based data. This means that a proper model for the galaxy surface brightness 
(SB) has to be convolved with the PSF and fitted to the observed SB. This is usually 
performed: 1) by one-dimensional fitting algorithms (e.g. Saglia et al. 1997) or 2) by 
two-dimensional methods (van Dokkum & Franx 1996; Simard et al. 1999). In case 1) the 
integrated light curve or the SB profile are derived by an elliptical fit to the galaxy isophotes 
and then fitted with the adopted model. In the two-dimensional approach, a direct fit of 
the galaxy SB is performed. Advantages and drawbacks of both methods have been largely 
discussed in previous works (see Kelson et al. 2000b, and references therein). In the present 
study we adopt the two-dimensional fitting approach for two reasons: 1) many galaxies in 
AC 118 belong to overlapping systems, and only the 2D approach can deal properly with 
these cases; 2) isophotal fitting of the ground-based images should be performed on a very 
small number of pixels, resulting in poorly constrained isophotes. 

The SB of distant galaxies is well described by two-dimensional models with constant 
cllipticity and position angle. In this study, we use two kinds of models: the Sersic model 
(Sersic 1968) and the (r 1//4 —) bulge + disk (B + D ) model. Since the B + D decomposition 
requires high resolution and signal to noise ratios (see Ratnatunga, Griffiths, & Ostrander 
1999), this model is applied only to the HST data. Ground-based structural parameters are 
derived by fitting the Sersic models. The models used for the PSF of the NTT and HST 
images are presented in § 3.1 while the fitting algorithm and the galaxy models are described 
in § 3.2. The limits of applicability of the fitting algorithm are set in § 3.3. 
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3.1. The Point Spread Function 

To model the PSF of the HST image, we used the Tiny Tim 5.0 software package (Krist 
1999). Since the PSF depends on the position on the field, different models were created 
for each galaxy at the corresponding chip location. Effects of PSF under-sampling and 
broadening from image jittering are smaller than few percents of r e for the typical sizes of 
galaxies at this redshift and therefore these were ignored (see Kelson et al. 2000b). 

To model the PSF of the NTT images we proceeded as follows. Stars in the cluster held 
were selected by the stellar index parameter SG of S-Extractor (Bertin & Arnouts 1996) and 
from the aperture magnitude. To minimize star - galaxy classification errors, we considered 
only objects with SG > 0.99. We selected sources with K < 18 mag in the SOF1 held and 
R < 21 mag in the region not covered from the infrared imaging. Such criteria allowed to 
obtain a sufficient S/N to model the PSF. The reliability of the selection is also supported 
by the photometric redshift algorithm: the colors of the candidate stars do not match galaxy 
template spectra at any redshift. 

The PSF of the NTT images is well described by a sum of two Gaussians convolved with 
the pixel function (see appendix B for details). The model has nine relevant parameters: the 
center coordinates, the central intensity, the Full Width at Half Maximum (FWHM) of the 
two Gaussians in the two perpendicular directions, the ratio of the central intensities, and 
the local value of the background. The fits were performed by a y 2 minimizing algorithm 
based on the Levenberg-Marquardt method (Press et al. 1986). Typically, the residuals from 
the fitting are smaller than 5%. The uncertainties on fitting parameters were estimated by 
numerical simulations on models with Gaussian noise. For each pixel, the standard deviation 
of the noise was computed by the distribution of the fit residuals in the adjacent pixels 2 . 
Figure 1 shows the values of the PSF parameters in the V —, R— and K —band as a function 
of the distance to the image centers. In general the PSF tends to become broader as this 
distance increases. The effect is very weak in the K —band image, while significant variations 
are found in the optical images. The trends of the parameters were modeled with empirical 
analytic curves, that were then used to obtain the PSF model relative to each galaxy position. 


2 This allows to take into account empirically the various sources of uncertainties: flat-field accuracy, 
photon and read-out noise, deviations from the fit model, etc. 



3.2. Modeling of the brightness distribution 


Galaxy images were fitted by the following model: 

M(BG, {p k }) — BG + B({p k }) o S , (1) 

where B is the galaxy brightness distribution, that is described by a set of parameters {pk}, 
S is the PSF model, BG is the value of the local background, and the symbol o denotes the 
convolution. For the Sersic models, the parameters {pk} are: the effective major semi-axis 
R e , the central surface brightness Jo, the Sersic index n, the axial ratio b/a , the position 
angle PA, the coordinates of the photometric center and the local value of the background. 
In the case of the B + D decomposition, the parameters that enter the fit are those of 
two Sersic models with n — 4 and n — 1 for the bulge and the disk respectively. The fit 
parameters are used to derive the total magnitude uit (e.g. Saglia et al. 2000), r e and 
{p) e - For Sersic models, r e is computed as the equivalent radius of the ellipse that encloses 
half of the total galaxy light, i.e. r e = \Jbfa ■ R e . For the B + D decomposition, r e is 
calculated from the growth curve of circular aperture magnitudes. In both cases, we have 
(/i) e = 2.51og(27r) + 51og(r e ) + %. For the Sersic models we found that both methods of 
computing r e agree within 2%. Therefore, the difference of the two definitions can be ignored 
in the following analysis. 

The parameters {pk} and BG were derived by a least squares fit, by minimizing the 
function 

X 2 (BG,{p k }) = f({p k })- £(/,-Mi(BG,{ Pt })) 2 , (2) 

- 3 

where Ij is the observed brightness distribution, j is an index running on image pixels and 
/ is a penalty function introduced to avoid ‘not allowed’ regions of the parameter space 
(i.e. unrealistic values that can be found on local minima of y 2 ). The y 2 minimization 
was performed by the Levenberg-Marquardt algorithm. The convolution with the PSF was 
computed straightly in the real domain in order to avoid high frequencies problems when 
convolving galaxy models by the Fast Fourier Transform. Bad pixels, degraded regions and 
pixels contaminated by background objects, were excluded from summation in equation 2. 
To this aim, a masking image was created interactively for each galaxy. Overlapping objects 
were fitted simultaneously. 

Uncertainties on structural parameters are usually estimated either from the formal 
errors on fitting results or by numerical simulations. The formal uncertainties are derived 
by the topology of the y 2 function near to the solution and therefore they do not take into 
account the various sources of errors on the estimated parameters (e.g. Smith et al. 1997). 
Numerical simulations give more reliable uncertainties but they depend on the adopted 
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galaxy model and are much more expensive in terms of computation time. As a compromise 
between accuracy and amount of computation time, we used a first order approximation: 
the uncertainties on the fit parameters were estimated from the derivatives of the y 2 func¬ 
tion taking into account the noise in the galaxy image and the uncertainties on the PSF 
parameters. Details are given in Appendix C. 


3.3. Simulations 


To test the performance of the fitting procedure and to set the limits to its applicability, 
we created two sets of simulations based on Sersic r 1 /" models convolved with the NTT 
A—band PSF and with the typical noise of ground-based data. To take into account the 
uncertainties on the PSF, all the simulations were generated with the same PSF but the fits 
were performed with different PSF models whose parameters were generated in agreement 
with the corresponding uncertainties (see § 3.1). 

The first set of simulations consist in r 1//4 galaxies spanning the full magnitude range 
in the R —band (19 mag < uit < 23 mag). Two values of r e were considered: r e = 0.3" and 
r e = 0.6", that correspond to the sizes of a small and a typical galaxy in a cluster at z ~ 0.3 
(e.g. Kelson et al. 2000a). In Figure 2 we give the results of the simulations in terms of 
differences between the input model parameters and the output from the fitting procedure. 
The quantities shown are r e , the axis ratio b/a (both given as logarithms), %, (/x) e , n, and 
the ‘fundamental plane parameter’, FP = log r e — 0.3 x that is the combination of the 
structural parameters entering the fundamental plane. It can be seen that the scatter of the 
fit parameters around the ‘true’ values is relatively small for rrir <22 where it amounts at 
most to 25% in r e , 0.5 mag in (fi) e , 0.15 mag in uit, 15% in FP and 10% in b/a. The scatter 
of n is higher, and varies from An ~ 1 at rriT = 19 mag up to An ~ 3 at rriT = 22 mag. 
Objects fainter than rriT = 22 mag are affected by larger and systematic differences between 
input and computed parameters. 


The second set of simulations was created in order to test the fitting algorithm with 
respect to the galaxy size. To this aim, a wide range of values of r e was explored (r e = 0.05" to 
r e = 3") for three Sersic models with n = 1,4, 7 respectively. For each value of r e , a suitable 
value of (/x) e was derived from the Kormendy relation, whose slope and the zero point 
were calibrated a posteriori from the present sample. In Figure 3, we show the comparison 
between true and estimated parameters. An inspection of Figure 3 shows that the parameters 
are recovered with a relatively small scatter (see above) and without systematic errors for 
r e > 0.15" (logr e > —0.8) for the whole range of n values. The shape of light profile is 
subject to a small systematic error (An ~ 1) for n = 7 below r e ~ 0.5" (logr e ~ —0.3). We 


rs_/ 
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also notice that the scatter increases for large n values, clue to the effect of the uncertainties 
on the PSF for galaxies with more concentrated profiles. 

The above results were used to derive the selection of the samples in this study. All 
the samples include the galaxies with r e > 0.15" in each band; in the R —band only galaxies 
with R <21 mag were included as a conservative limit; in the V —band a more restrictive 
selection in magnitude (V <20.5 mag) was necessary due to the lower S/N, while no further 
restriction was needed for the K —band photometry. 


4. HST structural parameters 

Structural parameters were derived for N = 50 galaxies in the HST image (sample A 
of Table 2), 42 of these objects are cluster members according to the photometric redshifts. 
Galaxies were fitted by both Sersic and B + D models, and a flag F was assigned to each 
galaxy for the different cases that were encountered during the fitting. For N = 5 galaxies, 
both models gave small fitting residuals, but the B + D decomposition produced meaningless 
values of the structural parameters for the bulge or for the disk component (e.g. r e < 0.02"). 
For such galaxies a possible second physical component is too weak or too small to be 
detected. However, except for one object that was excluded from the analysis, acceptable 
global parameters (r e and (/x) e ) were obtained from the B + D fit so that such galaxies were 
retained in the sample and were flagged with F = —1. For N = 24 objects, both Sersic 
and B + D models gave small fitting residuals and acceptable estimates of the parameters. 
They were flagged with F — 0. F = 1 was used to indicate N = 14 galaxies for which 
residuals from the Sersic model shown evidence for a second component that was properly 
described with the B + D decomposition. The flags F — 0 and F — 1 were assigned by visual 
inspection of the galaxy fits as illustrated in Figure 4. Finally, N = 7 galaxies, of which some 
objects with very distorted morphology and some BCGs with very large residuals (> 50%), 
were flagged with F = —2 and excluded from the present analysis. Fitting models suitable 
for such objects will be discussed in a forthcoming work. 

In Figure 5, we compare the structural parameters r e , (/x) e , rriT and FP obtained by 
different models. Table 3 shows the medians and the standard deviations of the differences 
in fitting parameters, and the relative bootstrap confidence intervals. By weighting each 
point in Figure 5 with the fitting uncertainties (§ 3.2), we did not find significant variations 
of the estimates given in Table 3. This also holds for the other comparisons of structural 
parameters in the present work (see § 4 and § 5). 

It is clear from Table 3 and from Figure 5 that the mean differences of the parameters 
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have no significant biases, also when galaxies with different flag F are considered. On the 
other side, the scatter in the structural parameters depends on the F value: as expected, 
minimum values are obtained for F = 0 when B + D and Sersic models give an equivalent 
description of the galaxy image. The scatter varies from 35% to 50% for r e , from ~ 0.8 
mag to 1.3 mag for (fj) e , and from ~ 0.25 mag to ~ 0.43 mag for my. One important 
fact is the robustness of the FP parameter against the adopted fitting method: due to the 
correlation between the errors in r e and in (/x) e (Hamabe & Kormendy 1987; Jprgensen, 
Franx,& Kjaergaard 1995), its scatter varies from ~ 7% for F = 0 to ~ 14% for F = —1. 
The previous results are fully consistent with the findings of Kelson et al. (2000a), although 
they find a smaller dispersion (~ 3%) between the FP parameters derived by different 
fitting methods. The results shown in Figure 5 and Table 3 give a measure of the intrinsic 
uncertainties on the structural parameters of the galaxies in the present sample and will be 
used as a reference in the following sections. 


5. Comparison of HST and NTT surface photometry 

The 2D fits for a representative subsample of galaxies in the HST and NTT R— band 
images are shown in Figure 6. Although the residuals from the fit are more evident in the 
HST images due to the higher resolution and to the higher S/N, both HST and ground- 
based images can be suitably fitted with the r 1 !' n models, also when overlapping objects are 
considered. 

To test the reliability of the ground-based structural parameters, we compared the 
estimates obtained by fitting Sersic models to the HST and the NTT images of galaxies in 
the sample B of Table 2. Since both filters cover approximately the same spectral range 
(rest-frame V— band), we can rule out any difference between the two datasets due to a 
possible waveband dependence of the structural parameters. The structural parameters are 
compared in Figure 7, where we plot the same quantities as in Figure 5. The medians and 
the standard deviations of the differences in logr e , (ju) e , m,T and FP are given in Table 4. 
Table 4 and Figure 7 prove that the mean differences of the parameters have no significant 
biases. The scatter between the two sets of measurements is relatively small: it amounts to 
~ 38% for r e , ~ 0.7 mag for ~ 0.3 mag for uit and ~ 13% for the FP parameter. 
It is worth to stress again that the FP parameter has significantly smaller dispersion with 
respect to the other quantities. The comparison with the results of the previous sections 
proves that the differences in the parameters derived from HST and from NTT photometry 
are of the same order of those found when applying different fitting methods to the HST 
image. In other words, the difference between the structural parameters derived from HST 
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and from NTT data is consistent with the precision with which these quantities can be 
measured (see Kelson et al. 2000a; Jprgensen et al. 1999). We thus conclude that reliable 
structural parameters can be still derived for typical sizes and luminosities of cluster galaxies 
at redshift 2 ~ 0.3 by using ground-based data taken in ordinary observing conditions (i.e. 
seeing ~ 1" and pixel scale r-u Q # 28"/pxl). 


6. Optical and Near-Infrared Structural Parameters 
6.1. Structure of Galaxies at z ~ 0.3 

The V — and I?—band structural parameters for N = 58 galaxies (sample C of Table 2) 
are compared in Figure 8 and Table 5. To allow a direct comparison, (/i) e and rriT of V— band 
galaxies were corrected with the observed V—R colors (Busarello et al. 2002b, in preparation). 
V — and R —band parameters show no significant systematic differences and the scatter 
between the two sets of data is of the same order as those found in the preceding sections, 
i.e. equal to the typical uncertainty of structural parameter measurements. The mean ratio 
t e,R,/re,v ~ O-^Tagi is consistent with that found in local early-type galaxies (r e j/r e y ~ 0.91, 
Scodeggio 2001; see also Scodeggio et al. 1998). We defer further considerations to § 7, where 
the results will be discussed in terms of color gradients. Our R— band structural parameters 
can be directly compared to those derived at z = 0.33 by Kelson et al. (2000b) in the same 
wavelength range. The full consistency of the two data sets is apparent from Figure 9 (upper 
panel) where the distributions of half-light radii are compared (the surface brightness is not 
shown for brevity reasons). 

The R— and IF—band structural parameters for the 93 galaxies in sample D (Table 2) 
are compared in Figure 10 and Table 6. NIR r e and (fx) e are significantly smaller than the 
relative optical parameters, while the n values are systematically higher. On the average, 
such differences amount to -0.125 in logr e , to -0.5 mag in (/j) e and to +0.5 in n. A deeper 
inspection reveals the existence of two distinct subsamples which are separated by the value 
of the Sersic parameter. For low n values (n < 4), the differences of K— to R— band effective 
radii amount to -0.09 in logr e , and the differences in the Sersic index are higher (An ~ 1 
instead of 0.5). On the other side, for n > 4 the Sersic indices are nearly the same in the two 
wavebands (An ~ 0) while the IF—band half-light radii are much smaller than in the optical 
(-0.2 dex). The confidence limits in Table 6 show that the separation into the two classes is 
highly significant. According to the Kolmogorov-Smirnov test, the probability that the two 
samples come from the same distribution is less than 2% for both An and Alogr e . 

The results of the previous sections make us confident that the systematic differences 
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found between the structural parameters are real and do not originate from biases induced by 
the fitting method. No other measurements of K —band (rest-frame H) structural parameters 
are available at z ~ 0.3 so that no direct comparison is possible. Again we defer further 
considerations to 5 7. 


6.2. Comparison with Local Galaxies 

To make a comparison with the structural properties of nearby galaxies, we made use 
of the sample of 70 galaxies belonging of the Coma Cluster from the compilation by Palire 
(1999). The distribution of the V —band r e from Pahre (1999) is shown in Figure 9 (notice 
that our R— band corresponds to the rest-frame V— band). The radii of Pahre (1999) were 
scaled to z — 0.31 by adopting Om = 0.31, Ha = 0 (different cosmologies do not significantly 
alter the result). 

The distributions of NIR r e are also compared in Figure 9 (the radii of Pahre 1999 were 
scaled as above). It is worth noticing that, although our rest-frarne waveband is different 
from the K —band in Pahre (1999), the two wavebands actually probe the same stellar 
populations, so that the direct comparison has an immediate physical interpretation. The 
medians of the distributions in the NIR are — 0.22 Iq 25 and —0.22Zq'\q for Coma and AC 118 
respectively, while the corresponding values in the optical are — 0.12 Iq 15 and — 0.08 Iq;° 2 - We 
thus conclude that the median of the half-light radii at z ~ 0.3 is fully consistent with the 
median at z ~ 0 both in the optical and in the NIR. 

This result indicates that the sizes of bright galaxies did not substantially change since 
z ~ 0.3 to the present epoch. Since the NIR probes directly the mass distribution inside 
galaxies, and assuming that changes in mass produce corresponding changes in size, we argue 
that during the last 4.4 Gyr 3 the bright cluster galaxies did not accrete a significant mass 
fraction and that (dissipation-less) merging did not play a major role (see also Nelson et al. 
2001 ). 


7. Implications for the Optical and Near-Infrared Color Gradients 

The role played by the Sersic index in the relation between NIR and optical structural 
properties (§ 6.1) suggests to discuss the combined behavior of half-light radii and Sersic 
indices in terms of internal color gradients. 


3 P m = 0.3, fi A = 0, H 0 = 50 km s^Mpc" 1 
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The available measurements of optical-NIR color gradients at intermediate redshifts 
concern 6 field early-type galaxies studied in V — H by Hinklcy & Im (2001). Unfortunately, 
no one of the galaxies is close enough in z to our sample to allow a significant comparison. The 
comparison of recent measurements of optical color gradients at intermediate-redshifts with 
those in nearby galaxies shows little or no evolution since z ~ 1, supporting the interpretation 
in terms of metallicity gradients (Peletier et al. 1990b,a; Saglia et al. 2000; Tamura & Ohta 
2000 ). 

We define the color gradient as A(/i 2 — /zi)/A(logr), where ^ is the surface brightness 
in the band i, and the differences are computed between an inner (r m ) and an outer ( tm ) 
radius. If the galaxy surface brightness in the waveband j is described by the Sersic law 

A ij(r) = R oj + 1-086 • bj{r/r e:j ) 1/nj , (3) 

we obtain 

A(/j 2 - Ui)/A(logr) = 2.5 log(e) • [log(r M ) - log(r m )] _1 • 

• [h ■ {{r m /r e , i) 1/ni - {r M /r e , i) 1/ni ) + b 2 • (( r M /r ej2 ) 1/n2 - (r m /r e , 2 ) 1/n2 )] , 

where b is a function of n (b ~ 2n — 1/3, see Caon, Capaccioli, & D’Onofrio 1993) and 
the indices 1,2 denote the two wavebands. In order to make a direct comparison with 
previous works (Peletier et al. 1990b; Saglia et al. 2000) we set r m = 0.1 - r e R and r M = r e R . 
Equation 4 is more suitable for the present analysis than the approximation of Sparks & 
Jprgensen (1993), in that it allows to address the case of different values of n and it is exact 
for arbitrary ratios of the half-light radii. We also verified that this formula gives a very 
close approximation to the linear fit to the color profiles. 

The V—, R— and K —band structural parameters derived in the previous sections were 
substituted in equation 4 to derive the implied color gradients. The uncertainties on the 
gradients were derived from propagation of the errors on r e and n, taking into account the 
covariance term. The V — R and R — K color gradients derived by means of equation 4 are 
shown in Figure 11 as a function of n, and the relative median values are given in Table 7. 
The average optical (V — R) color gradient is *-0.06l2‘i5 mag/dex, in agreement with the 
V — I gradients derived by Tamura & Ohta (2000) at z =0.37 and with the optical (V, R, I ) 
gradients by Saglia et al. (2000) at z ~ 0.38 for cluster early-type galaxies (see Figure 6 in 
Saglia et al. (2000) and Tables 2 and 3 in Tamura & Ohta 2000). 

The V — R gradients present a weak trend with n: galaxies with n < 4 have larger 
negative gradients, with an average of -0.19 mag/dex, while for n > 4 the average gradient is 
consistent with zero or slightly positive (+0.07 mag/dex). The wide majority of the galaxies 
have negative R — K color gradients: the median value for the whole sample is — 0.49 Iq 5 5 
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mag/dex. The trend with n is in this case stronger: the average R—K color gradient is —0.73 
mag/dex for n < 4 and —0.35 mag/dex for n > 4. Furthermore, the results of § 6.1 indicate 
that the color gradients are mainly due to the difference in n for n < 4 and to the difference 
in r e for n > 4. The optical-NIR color gradients are thus related to the relevance of the disk 
component, as it is apparent from Figure 11: galaxies with lower bulge-to-disk ratio (n —> 1) 
have large R—K color gradients, while bulge-dominated galaxies have moderate gradients. 

In order to probe the origin of color gradients with the large wavelength baseline provided 
by the K— band, we adopted two different models in which color gradients are driven by 
(A) ‘pure age’ or (M) ‘pure metallicity’ gradients. We used synthetic colors from the code 
of Bruzual & Chariot (1993) to simulate early-type galaxy spectra. An exponential star 
formation rate with time scale of 1 Gyr was adopted. In model (A) we assumed that the 
inner stellar population (r ~ r m ) is 6.25 Gyr old at z = 0.31. This value corresponds to the 
mean luminosity-weighted age of stellar populations in AC 118 galaxies as inferred from the 
color - magnitude relation (see Merluzzi et al. 2001). The stellar population age at the outer 
radius (r ~ tm) was adjusted to match the observed color gradients, and the metallicity was 
fixed to Z = Z Q . In model (M) the inner stellar population is characterized by a metallicity 
Z = 1.8 Z Q (Tamura et al. 2000) and the observed color gradients were used to derive the 
outer Z value. Age was fixed to 6.25 Gyr. From the median values of the R — K color 
gradients we derived the outer age and metallicity (tm and Zm) for models (A) and (M) 
respectively. Such values were then used to predict the V — R color gradients at z = 0.31 
and the V — R and R—K gradients at z ~ 0. 

Table 8 gives the predicted V—R and R—K gradients along with the gradients observed 
both at z ~ 0.3 and at z ~ 0. Optical-NIR color gradients were measured for 12 nearby 
bright elliptical galaxies by Peletier et al. (1990a, b). From their data we obtain an average 
R — K gradient of —0.14, with a broad distribution of values ranging from —0.62 to +0.15. 
Palire et al. (1998) estimated an average V — K gradient of -0.18 by applying the formula 
of Sparks & Jprgensen (1993) to the ratio of the relative half-light radii, while an average 

V — H gradient of ~ —0.18 can be inferred in the same way from the Coma cluster H and 

V effective radii given by Scodeggio et al. (1998). 

The results shown in Table 8 may be summarized as follows. The pure age model is 
able to explain the evolution of the R — K color gradient between z ~ 0.3 and z ~ 0, but 
it fails to describe the V — R gradient at z ~ 0.3. Conversely, the pure metallicity model is 
able to reconcile the V — R and R—K color gradients at z ~ 0.3, but it fails to describe the 
evolution of the R—K gradient as a function of look-back time. 
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8. Discussion and Conclusions 

In this study we derived optical and NIR structural parameters for galaxies belonging 
to the cluster AC 118 at z = 0.31 from ground-based photometry. The present data set con¬ 
stitutes the first large sample of structural parameters determined at intermediate redshifts 
in different wavebands and, in particular, it is the first such sample in the near-infrared. The 
samples consist in N=93 objects studied in the R and K wavebands and in N=58 galaxies 
studied in the V— band. 

The reliability of the data was verified by comparison with the results from HST pho¬ 
tometry for a subsample of galaxies: the scatter between measurements from HST- and from 
ground-based photometry is equivalent to the scatter between parameters determined with 
two different fitting methods from HST images, namely of ~35% in r e . Given that these 
values of the scatter are generally claimed to be a measure of the ‘intrinsic uncertainty’ in 
structural parameter measurements, we can state that reliable structural parameters can be 
still derived from ground-based photometry in ordinary observing conditions at least up to 
z ~ 0.3. 

Since morphology was not a selection criterion in this study, some remarks are needed 
for the comparison of our results with other works. Couch et al. (1998) established that 78% 
of the galaxies in the core of AC 118 are early-types while 22% are spirals. They also found 
a significant rise in the fraction of SO galaxies with respect to the local density-morphology 
relation. Morphological classification from Couch et al. (1998) is available for 21 galaxies 
in sample D. The number of spirals is 5, of which 3 have n <1.5 while the other two have 
n > 4. Only one of the spirals has ‘late-type’ colors, while the others lie on the red sequence 
of AC 118. It is interesting to notice that all the remaining galaxies with ‘late-type’ colors 
lie outside of the cluster core and have n < 2. They are likely disk galaxies dominated by 
young stellar populations. These facts are in agreement with a scenario in which spirals are 
accreted into the cluster from the held, cease to form stars, become gas-deficient objects 
and eventually undergo a morphological transformation (Couch et al. 1998). The remaining 
galaxies with morphological information consist in 10 SO, eight of which with 2 < n < 3 and 
two with n ~ 8, and 6 ellipticals, for which we find n > 4. We stress that the galaxies with 
large disks and late-type colors are only 4% of galaxies in sample D. 

The main feature concerning the multi-band properties of galaxies is that these are much 
more concentrated in the NIR than in the optical. The mean ratio of NIR to optical half-light 
radii (~ 0.75) is not different from that found in local samples of early-type galaxies (~ 0.8, 
e.g. Pahre 1999), but the situation is complicated by the systematic differences found in 
the shape of the brightness distribution: galaxies with low n value have a light profile more 
peaked in the center in the /7-band (n K — Ur ~ 1), while for n > 4 a lower r eK /r ej R is 
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observed (~ 0.6). Such trends can be interpreted as due to the increase of the disk fraction 
when n approaches unity. 

To discuss the combined effects of half-light radii and Sersic indices, we derived the 
optical and NIR color gradients by means of equation 4. Most of galaxies in AC 118 do not 
show significant optical color gradients, with a median value of — 0 . 061 q'i 5 - This result is 
consistent with previous estimates of optical color gradients at intermediate redshifts (Saglia 
et al. 2000; Tamura & Ohta 2000). Conversely, strong negative gradients are found in optical- 
NIR colors: the median of the R—K gradients amounts to —0.49 Iq 55 . Moreover, this effect 
depends on the shape of the light profile: the R — K gradient increases from — 0 . 73 log 3 for 
galaxies with low n value to —0.35 Iq 4 ® for n > 4. Again this result can be explained by an 
increase in the disk fraction when n —> 1 . 

To interpret these results in terms of the properties of the stellar populations, we intro¬ 
duced a pure age and a pure metallicity model. We used the R — K gradients to calibrate 
both models and to predict V — R color gradients at z = 0.31. Moreover, the evolution of 
color gradients as inferred by the models was compared with the measurements by Peletier 
et al. (1990a,b) for a sample of nearby elliptical galaxies. The age model is not able to 
reconcile our optical and NIR estimates at z — 0.31 since a significant negative gradient is 
predicted in V — R (~ —0.26). On the other hand, the predicted color gradient agree with 
the measurements at z = 0. Conversely, the pure metallicity model gives correct results at 
z = 0.31 but predicts too strong NIR gradients at z = 0. Since the age model is able to 
describe the evolution of color gradients, whereas the metallicity model correctly explains 
the dependence of color gradients on wavelength, we argue that a model in which effects of 
age and metallicity are consistently combined could explain optical and NIR color gradients 
from z = 0.31 up to z = 0. Finally, we remark that dust likely affects the optical-NIR color 
gradients. Dust is actually presents in the core of most early-type galaxies (Goudfrooji & de 
Jong 1995; van Dokkum & Franx 1996) and is known to produce high central color gradients 
in bulges of SO and spiral galaxies (e.g. Peletier et al. 1999). 

The distributions of both the optical and NIR half-light radii of galaxies at z ~ 0.3 
are found to be consistent with those of Coma early-type galaxies. The result for the NIR 
waveband is the most significant, since the NIR probes directly the mass distribution in¬ 
side galaxies. The observed distributions thus indicate that the sizes of bright galaxies did 
not substantially change since z ~ 0.3 to the present epoch. This result suggests that 
(dissipation-less) merging did not play a major role in galaxy evolution during the last 4.4 
Gyr (see also Nelson et al. 2001). 


We are grateful to Adam Stanford who provided us with his own photometry of AC 118, 
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helpful comments and suggestions and Dr. M. D’Onofrio helpful discussions. The observa¬ 
tions at ESO were collected during the guaranteed time of the Osservatorio Astronomico di 
Capodimonte. This work is partially based on observations made with the NASA/ESA Hub¬ 
ble Space Telescope, obtained from the data archive at the Space Telescope Science Institute. 
STScI is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract NAS 5-26555. Michele Massarotti is partly supported by a MURST-COFIN 
grant. 
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A. Calibration of HST data 

AC 118 was observed on November 1994 with the HST Wide-Field Planetary Camera 
2 in the F702W filter. Two exposures of 2300 s and one exposure of 1900 s for a total 
integration time of 6500 s are available for a single pointing. 

The HST images were calibrated into the Johnson R —band by means of the NTT 
photometry. We proceeded as follows. The frames were smoothed to the same resolution of 
ground-based images by a convolution with the NTT R —band PSF. Aperture photometry 
was then performed with the SExtractor software (Bertin & Arnouts 1996) in the same 
aperture of the ground data. Figure 12 shows the difference of R and F702W instrumental 
magnitudes versus V — R colors. The data points were fitted with the equation: 

R = F702W instr A a * (V — R) + b . 

We found a = 30.93T0.10 and b = 0.30A0.12, and that the root mean square of the residuals 
to the fitted relation varies from 0.05 mag for R < 20.5 mag to 0.1 mag for the whole sample. 
The few galaxies (9%) which deviate more than ~ 0.2 mag from the fit, turned out to be 
objects contaminated by bright neighbors. In the same figure we also plot the synthetic 
transformation to the R waveband of Holtzman et al. (1995, see their equation 9 and Table 
10 ): 


R = F702W instr A 0.486 * (V - R) — 0.079 * (V - R) 2 i 

i 2.5 * log (t exp ) i 2.5 * log (GR) A A m ap A A m exp ? (Al) 

where t exp is the exposure time, GR is the gain ratio of the WFPC2 chips, A m ap is an 
aperture correction term derived from the encircled energy measurements of Holtzman et 
al. (1995), and A m exp is the long exposure correction term from Hill et al. (1998). The 
same gain ratio GR = 2 was adopted for each chip, since the effects of gain differences are 
very small (< 0.01 mag) when compared to the other sources of uncertainty. We adopted 
A m ap = AO.10 mag and A m exp = A0.05 mag. A further shift of A0.05 mag was needed 
to better match equation Al to our data. This is in agreement with previous results of 
calibration of HST images with ground-based photometry (see Jprgensen et al. 1999). 


B. PSF of ground-based images 


The PSF of NTT ground-based images was described by the following model: 

S(x,y) = [Gi + G 2 ]*R(x,y) , 


(Bl) 
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where S is the PSF brightness distribution, Gj{j = 1,2) are two-dimensional Gaussian func¬ 
tions and R is a box function equal to 1 on the pixel area and 0 elsewhere. The convolution 
in equation B1 was computed from the one-dimensional convolution of a Gaussian with a 
rectangular function. For the one-dimensional case, the following expression is obtained: 


g*r (x) — (2a) 1 



+ Erf 


a/2—x 
y/2 a 


(B2) 


where g(x) is is a normal Gaussian function of width parameter a, r{x ) is a rectangle of 
width a, and Erf(x) is the error function. The analytical expression of the ground-based 
PSF is obtained by substituting the formula B2 in equation Bl. 


C. Uncertainties on galaxy structural parameters 

Galaxy structural parameters were derived by minimizing the right side of equation 2, 
i.e. by solving the following system of equations: 

^G 2 ({ft} I {h }. M) = 0 . (Cl) 

where y 2 is a function of the fitting parameters {pk}, {Ij} are the brightness values of image 
pixels, and {s;} are the PSF parameters for the ground-based data. The changes {Slj}, 
of {Ij} and {s/} determine variations {8pk} of the solutions {pk} of equations Cl. At 
the first order, {8pk} are obtained from the following formulae: 

Sp = H~ l ■ (A ■ SI + B ■ 8s) , (C2) 

where the symbol • denotes matrix multiplication, 8p, 81 and 8s are column vectors of 
components 8pk, 8Ij and 8si respectively, H~ 1 is the inverse of the y 2 Hessian matrix, A 
and B are rectangular matrices of components <9y 2 /dpkdsi and <9y 2 /dpudlj respectively. All 
quantities are computed at {pk}, {Ij} and {sj}. 

The noise of image pixels and the uncertainties on PSF parameters (see § 3.1) were used 
to estimate the quantities {8Ij} and {5s;}. Uncertainties on fitting parameters were then 
derived by equations C2. 
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Fig. 1.— Dependence of the gronnd-based PSF on the chip position. The quantities shown 
in the figure are: the ratio of the central intensities of the ‘outer’ and ‘inner’ Gaussians / 2 // 1 , 
the FWHM in the x- and y- directions of the inner Gaussian (FWH and FWHMy l> ), 
the FWHM in the x- and y- directions of the outer Gaussian (FW H and FWHMy 2) ) } 
and the asymmetry A = FWHM^/FWHMy 1 ' 1 of the inner Gaussian. On the x-axis the 
distance to the frame centers d is plotted. For each quantity, three panels are shown for the 
V —, R— and K —bands as indicated in the npper-left plots. The error bars correspond to 
la standard intervals and the solid lines are the curves used to create the PSF models for 
each galaxy position (see text). 

Fig. 2.— Results of the numerical simulations for R-band photometry. The differences 
between input and output values are shown for the following quantities: logarithm of r e , 
(fj) e , itt-t, fundamental plane parameter FP, logarithm of the axis ratio b/a, and n. For each 
quantity, the two panels correspond to simulations with r e = 0.3" and r e = 0.6". The total 
magnitude of the models is given on the abscissae. 

Fig. 3.— Minimum r e for the applicability of the 2-D fitting to the ground-based data. The 
same quantities as in Figure 2 are plotted on the vertical axis while the logarithm of the ‘true’ 
effective radius is shown on the abscissae. For each quantity, the three panels correspond to 
the different values of n (see upper-left panels). 

Fig. 4.— Two-dimensional fit of galaxies in the HST image by means of Sersic and of B + D 
models. Left and right pictures are relative to galaxies with F = +1 and F — 0 (see text). 
Upper panels: original image; middle left and right panels: Sersic and B + D fitting models; 
lower panels: corresponding residuals. 

Fig. 5.— Distribution of the differences between structural parameters derived from the 
Sersic fit and from the B + D decomposition for HST galaxies. From left to right, and 
from top to bottom: r e , (/x) e , uit and FP parameter. Solid, dashed and dot-dashed lines 
are relative to galaxies with different fitting flag F as indicated in the upper-left panel. 
Differences are in the direction Sersic-( B + D). 

Fig. 6.— Two-dimensional fit of galaxies in the HST F702W and NTT R —band images. For 
each quadrant, the left and right panels correspond to NTT and HST data respectively, and, 
from top to bottom, galaxy image, fitting model, and image of residuals are shown. HST 
and NTT images are scaled to the same counts/arcsec 2 . 
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Fig. 7.— Comparison of structural parameters derived from NTT R and HST F702W images. 
For each panel, the histograms show the fractional distribution of differences between NTT 
and HST parameters. The plotted quantities and the plotting ranges are the same as in 
Figure 5. Differences are in the direction NTT - HST. 


Fig. 8.— Differences (V — R ) of structural parameters derived from V and R images. The 
R —band n is plotted on the abscissae. 


Fig. 9.— Distribution of optical (upper panel) and NIR (bottom) half-light radii of: 1) 
cluster galaxies at z = 0.31 (solid line, this work); 2) cluster galaxies at z = 0.33 (dashed 
line, Kelson et al. 2000b); early-type galaxies in Coma (z = 0.023, dash-dotted line, Pahre 
1999). The Coma radii were transformed to z — 0.31 by adopting Dm = 0.3 and Da = 0. 


Fig. 10.— Differences (K — R) of structural parameters derived from R and K images. The 
R —band n is plotted in the abscissae. The plotting ranges are the same as in Figure 8. 


Fig. 11.— Color gradients of AC 118 galaxies as derived from equation 4 (see text). V — R 
and R — K color gradients are shown in the upper and lower panel respectively. The Sersic 
index in the R —band is plotted on the abscissae. Error bars define lcr standard intervals. 


Fig. 12.— Calibration of HST data by ground-based NTT photometry. Differences of R and 
F702W instrumental magnitudes are plotted versus V—R. The solid line is the least squares 
fit to the data. The dashed curve represents the synthetic transformation of Holtzman et al. 
(1995). 
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Table 1. Summary of the photometric data. 


Band 

Instrument 

Exp. Time 
(sec) 

Field 

pixel scale 
(arcsec/pxl) 

seeing 

(arcsec) 

V 

NTT-EMMI 

2700 

9/15x8/6 

0.267 

1.2 

R 

55 

1800 

55 

55 

1.0 

K 

NTT-SOFI 

15900 

5'x5' 

0.292 

0.8 

F702W 

HST-WFPC2 

6500 

5.8 arcmin 2 

0.1 
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Table 2. Description of the samples. 



A 

B 

c 

D 

wavebands 

F702W 

R, F702W 

V, R 

ii, K 

phot, redshift 

all 

all 

0.3 

0.3 

selection 

R < 21, r e > 0.15" 

R < 21, r e > 0.15" 

R <21, V < 20.5, r e > 0.15" 

R < 21, r e > 0.15" 

N 

50 

37 

58 

93 

fitting models 

B + D, r 1 / 71 

r l/n 

r l/n 

^*1 j rr 






Table 3. 


Comparison of Sersic and B + D fitting parameters. 


F 

Alog 

r e 

A(/i) e 


Amy 


AFP 



m 

std 

m 

std 

m 

std 

m 

std 

-1 

- 0 . 06 t°;°l 

a oo0.47 
U,OZ 0.29 

n 40+ 0,24 
u.4u_ 0-90 

l.^i 2 

—0 If) - * -0 ' 26 
u>lu -0.23 

n aqO.95 
u -^' 3 0.37 

_i_0 06~*~°' 10 

-i-u.uo +0 00 

f) 1 fi 014 
u - iu 0.09 

0 

0.002l°;°24 

o.isgdi 

0 Q3+ 0 - 14 

u -uo_ 0 .io 

n cft0.68 

u,oo 0.50 

0 o4+ 0 - 07 
U - U4 _0.001 

n i 7O.21 
u - x '0.15 

-0.0ll±g;ggi 

f) 047 0 - 057 
U - U4: ' 0.042 

1 

-o. 06 ± 8;“ 3 

0.208:?! 

1 

o 

to 

o 

1 + 
o © 

^ o 

O -SI 

n 7Q1-20 
u - ‘ y 0.66 

fi -|2+ 0 - 18 
u - lz +0.04 

n 04 0 34 
u - Z4: 0.21 

— 0 . 003 ±ggg 5 

0 f)3fi 0 060 
U.UoOq 030 

-1,0,1 

—0 n^+ 0 - 001 
u.uO—0 056 

0 - 1858 :?? 

— 0 1 qTO.016 
U - 1U -0.21 

0 7fi 0-87 
u - '°0.70 

0 04” 1 " 0 ' 08 
U - U4 +0.006 

0 24 0-28 
u - z m.22 

O.OOljlg.go? 

n n«0.0 7 
u - ud 0.05 


Note. — Column 1: fitting flag F. Columns 2, 3, 4 and 5: medians and standard deviations of the differences of the structural 
parameters. The upper and lower values define the 68% confidence interval. Differences are in the direction Sersic - (B + D). 
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Table 4. Comparison of HST and NTT R —band fitting parameters. 



m 

std 

A(log r e ) 

( N 

O 00 

o o 

6 o 
+ 1 

o 

o 

1 

0 1 fi 0-21 
U - 1D 0.14 

a ((mM 

— fi 1 004 

u ’ lo -0.26 

oo 

000 

do 

t- 

cq 

o 

A(mr) 

CO IN 

O tH 

d o 

1 1 

00 

o 

o 

1 

o so 0-33 
u - t)u 0.28 

A (FP) 

-i-O 02f)”*~ 0 ' 037 
^u.uzu +0 004 

0 Ofi 0 08 
u>ud 0.05 

A (log b/a) 

+0-Oll|aoo3 

0 082 0 095 
U.UoZq 075 

A(n) 

+o-oitS;Jo 

2 fi 2,9 
z,0 2.4 


Note. — Median and standard devia¬ 
tion of quantities in column 1 are shown in 
columns 2 and 3. Upper and lower values 
define the 68% confidence interval. Differ¬ 
ences are in the direction NTT-HST. 
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Table 5. Comparison of NTT V— and R —band structural parameters. 


m std 


A (log r e ) 

+ 

© 

o 

4^ 

1 + 
o © 

o o 
o 

4+ 

0 17O.26 

A ((+)e) 

+0.019i°;J 2 7 

f) 7KI.OO 
u - '°0.63 

A (raj>) 

-0.045I°;°83 

0 26 0 34 
u - zo 0.22 

A (FP) 

+0.008to.oo? 

O 

O 

00 

00 

bo 

-4© 

A (log b/a) 

-0.0101°;°?? 

r-m 

00 

do 

CO 

O 

O 

A(n) 

+0.0241°;?? 

1 Q 2 - 2 


Note. — See Table 4. Differences are 
in the direction V — R. 







Table 6. 


Comparison of NTT K— and R —band structural parameters. 



all 


n < 4 


n > 4 



m 

std 

m 

std 

m 

std 

A(logr e ) 

-0.1251°;“ 

o 1 qO- 21 
u - iy 0.18 

-o.o9:°;?I 

n 12°- 14 

U.1Z 0 11 

—0 20 _0 17 
u.zu_ 0 24 

0 2S°- 26 

u.zoq.21 

A((M>e) 

-omzUI 

n q4l.02 
u,y ^0.89 

-0.371“;“ 

0.68g;|5 

fl 7Q — 0-66 

U * '°-0.93 

1+2 l;i 

A(mr) 

+0.04i°;°» 

n 2fi 0 29 
u - zd 0.25 

—0 04“ 9 02 
u.u4_ 0 .io 

o.isgfi 

_|_f) I k+O- 2 ^ 

-|-u.io_|_ 0 os 

0.29°;1? 

A (FP) 

_i_0 04+ 0 - 05 

0 11 012 
u - iJ -o.io 

+0-04+°;°?5 

0 f)Q 011 
u - uy 0.08 

_i_f) qoc+ 0.053 
i-u.uoo +0 015 

0 1 s 014 
u - lci 0.12 

A(log b/a) 

—0 002” 1 ” 0 005 
u.uuz_ 0 009 

0 OfiO 0,067 
U.UOOq 055 

-0.005t°°°I 

0 OfiO 0,067 
u.uoOq 056 

r- co 

0 0 

0 0 
d d 
+ 1 

0 

0 

O 

d> 

+ 

0 BO 0,76 

U - OU 0.52 

A(n) 

_|_f) co - *" 0-60 
-t-u.oz +0 48 

2 0 2 - 2 
z,u 1.9 

_i_0 95"*" 1,0 
-|-u.yo_|_ 0 g 5 

1 6 1 - 9 

+0.08+°;“ 

9 92.5 

z.z 2 0 


Note. — Medians and standard deviations are shown for galaxies with different ranges of n. Differences are in the 
direction K — R. 
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Table 7. Median values of V — R and R—K color gradients of AC 118 galaxies. 




all 

n < 4 

n > 4 

(V- 

R) 

-o.oetn? 

—0 IQ -0 ' 09 
u - iy -0.28 

+0.07tg;J® 

(R- 

K) 

-0.49 = °;« 

—0 73 —0,63 
u - ‘°-0.83 

-0-35:°;” 


Note. 


Color gradients are in units of mag/dex. 







Table 8. Computed and observed color gradients. 


MODEL 

grad(V — R) 
z ~ 0.3 

grad(V — R) 
z ~ 0 

grad(R — K ) 
z ~ 0.3 

grad(R — K ) 

2 ~ 0 

AGE 

-0.26 

-0.04 

-0.49 a 

-0.13 

TM = 3.2 Gyr 

METALLICITY 

-0.15 

-0.10 

-0.49 a 

-0.41 

Z M = 0.5 Zq 

observed 

-0.1 b 

-0.02 c 

-0.49 d 

-0.14° 


a Constrain to the model. 

b See Tamura & Ohta (2000); Tamura et al. (2000); Saglia et al. (2000), and this 
work. 

c See Peletier et al. (1990a, b). 

d This work. 

Note. — The predictions are based on pure age and pure metallicity gradient 
models (see text). Column 1: adopted model. Column 2, 3: optical - optical color 
gradient at z ~ 0.3 and z ~ 0 respectively. Column 4: optical - NIR color gradient 
at z = 0. color gradients are reported in units of mag/dex. 
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